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We perform a new study of future sensitivities of galaxy redshift surveys to the free-streaming 
effect caused by neutrino masses, adding the information on cosmological parameters from measure- 
ments of primary anisotropies of the cosmic microwave background (CMB). Our reference cosmo- 
logical scenario has nine parameters and three different neutrino masses, with a hierarchy imposed 
by oscillation experiments. Within the present decade, the combination of the Sloan Digital Sky 
Survey (SDSS) and CMB data from the PLANCK experiment will have a 2a detection threshold on 
the total neutrino mass close to 0.2 eV. This estimate is robust against the inclusion of extra free pa- 
rameters in the reference cosmological model. On a longer term, the next generation of experiments 
may reach values of order ^ = 0.1 eV at 2cr, or better if a galaxy redshift survey significantly 
larger than SDSS is completed. We also discuss how the small changes on the free-streaming scales 
in the normal and inverted hierarchy schemes are translated into the expected errors from future 
cosmological data. 

PACS numbers: 14.60.Pq, 95.35. +d, 98.80.Es 



I. INTRODUCTION 

Neutrino physies has provided the first clear indication of particle physics beyond the Standard Model, since we 
have experimental evidences for non-zero neutrino masses. Analyses of data from atmospheric and solar neutrino 
experiments have shown the allowed regions for the squared mass differences (Am^) at two different scales. Such 
values will be known with better precision in the next years, in particular for the larger atmospheric AmJ using the 
results of future long-baseline oscillation experiments. 

However, from oscillation experiments no information can be obtained on the absolute values of neutrino masses, 
since the lightest neutrino mass remains unconstrained. Tritium decay experiments tell us that each neutrino mass 
cannot be larger than 2.2 eV (95% CL) at present [H, to be improved to ~0.35 eV with KATRIN 0. More stringent 
bounds exist from experiments searching for neutrinoless double beta decay, that will be improved in the near future 
0, but unfortunately they depend on the details of the neutrino mixing matrix. 

Cosmology offers several advantages: the cosmic neutrino background provides an abundant density of relic neu- 
trinos with an equal momentum distribution for all flavors (up to 1% corrections), which implies that mixing angles 
have no effect. Although neutrinos cannot be the dominant dark matter component, they can still constitute a small, 
hot part of the matter density producing an erasure of perturbations at small scales through their free-streaming 
effect (for a review, see e.g. 0). A comparison with data from the large scale structure (LSS) of the Universe is thus 
sensitive to neutrino masses, as emphasized in Q. 

At present, cosmological data allow us to bound the total neutrino mass to values of ^rrii, < 0.6 — 1.0 eV 
d 0, H H Efl m 113 , depending on the data and priors used. These ranges already compromise the 4 neutrino 
scenarios that could explain the additional large neutrino mass difi'erence required by the LSND results (that also 
imply a fourth, sterile neutrino), but is not yet capable of reaching the necessary 0.1 eV range in order to test the 
hierarchical 3 neutrino schemes. But such small masses could be detected in the next future when more precise 
cosmological data are available, in a parallel effort to those of beta and double beta decay experiments on Earth. 

In this paper we analyze the future sensitivities of cosmological data to neutrino masses, extendin g th e pioneering 
work and in particular the detailed analysis in (see also Q)i that was more recently updated in |l5l|. In contrast 
to this last work we consider, in addition to ideal Cosmic Microwave Background (CMB) observations limited only 
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FIG. 1; The two neutrino schemes aUowed if Amatm ^ ^^aun- normal hierarchy (NH) and inverted hierarchy (IH). 



by cosmic variance, the experimental specifications of satellite missions such as PLANCK and the mission concept 
CMBpol (Inflation Probe), as well as ground-based detectors such as ACT or SPTpol, that will extend the PLANCK 
data to smaller angular scales. We also increase the number of cosmological parameters of previous analyses, including 
also the helium fraction, extra relativistic degrees of freedom, spatial curvature, dark energy with constant equation 
of state, or a primordial spectrum with running tilt. Finally, our work is the first one in which it is assumed that 
neutrinos have three different masses, in order to compute accurately the free-streaming effect associated to the mass 
schemes allowed by oscillation experiments. 

Note that throughout this work, we will assume that the LSS power spectrum is measured solely with galaxy 
redshift surveys. For complementary constraints based on gravitational lensing, we refer the reader to Refs. Il7j. 

This paper is organized as follows. In Sec. II we review the expected values of neutrino masses and their impact on 
Cosmology. We describe future CMB experiments and galaxy surveys in Sec. Ill and the method to forecast the errors 
on cosmological parameters in Sec. IV. Finally, we present our results in Sec. V, with a summary and conclusions in 
Sec. VI. 



II. NEUTRINO MASSES 



Nowadays we have experimental evidences for neutrino oscillations from solar and atmospheric neutrino detectors, 
recently also supported from data on neutrinos from artificial sources (Kamland and K2K). Detailed analyses of the 
experimental data lead to the following values of the mass squared differences (best fit values ± 3a ranges) 

Aml,^ = Ami, = (2-6tl:^) x IQ-^ eV^ 

AmL. = Ami, - (6.9t?J) x 10"^ eV^ (1) 

taken from These ranges are only slightly different in other recent analyses, see e.g. |23|, while a lower 

^'^atm seems required by new Super-Kamiokande data and 3-dimensional atmospheric fluxes. The errors in the 
above equation will be significantly reduced with new data from Kamland in the case of Am^i, and with data from 
future long-baseline oscillation experiments such as MINOS, ICARUS and OPERA, which will give the atmospheric 
Am^ with 10% accuracy (reduced to 5% with the superbeam proposal JPARC-SK) Current data also provide 

the allowed ranges of the neutrino mixing angles 612 and 023, and an upper bound on ^13. 

Indications for a third, heavier Ato^ exist from the LSND experiment p3 |. implying a fourth (sterile) neutrino. 
Such a mass is already being tested by present cosmological data, although not ruled out yet 0, S 0> ' ^^'^ the 
LSND results will be checked by the ongoing experiment MiniBoone. Here we choose not to include such a large Amf, 
and consider only the values in Eq. ^ 

The three neutrino masses that lead to the values in Eq.^can be accommodated in two different neutrino schemes, 
named normal (7713 > m2 > mi) and inverted (to2 > 'm-i > 'm-s) hierarchy, as shown in Fig. ^ that we will denote 
NH and IH. At present we have no indication of which scheme is the correct one. However, it has been suggested 
that some information could be extracted from future data from Supernova neutrinos, very large baseline oscillation 
experiments, or neutrinoless double beta decay searches if the effective mi, is below some threshold (for reviews, see 
6-g- |23 . l24j '). In general, determining the type of mass spectrum depends on the precision with which the other mixing 
parameters would be measured. 

Relic neutrinos were created in the Early Universe and decoupled from the rest of the plasma when the temperature 
dropped below ~ 1 MeV, when they were ultra-relativistic. After decoupling all neutrino flavors kept a Fermi-Dirac 
spectrum, only distorted at percent level during the process of electron-positron annihilations into photons psl l26l | . 
It is well-known that massive neutrinos could account for a significant fraction of the total energy density of the 
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FIG. 2; Neutrino masses as a function of the total mass in the two schemes for the best-fit values of Am^ in Eq. The 
vertical line marks the smallest value of ^ in the inverted scenario. 



Universe today, being their contribution directly proportional to the number density. For vanishing neutrino chemical 
potentials, the total neutrino contribution to the critical density is given by 

93.2 cV ' ^ ' 

where h is the Hubble constant in units of 100 km s^^ Mpc~^ and ^ rrii^ runs over all neutrino mass states. For fixed 
neutrino masses, fli, vifould be enhanced if neutrinos decoupled with a significant chemical potential (or equivalently, 
for large relic neutrino asymmetries), but this possibility is now ruled out |27l |. 

Therefore cosmology is at first order sensitive to the total neutrino mass ^ — mi + + ™3 (for the 3 
neutrino schemes that we consider), but blind to the neutrino mixing angles or possible CP violating phases. This 
fact differentiates cosmology from terrestrial experiments such as beta decay and neutrinoless double beta decay, 
which are sensitive to \Uei\'^m^ and | U^^mi\, respectively, where U is the 3x3 mixing matrix that relates the 
weak and mass bases. 

It is interesting to see how the total mass is distributed among the neutrino states for the two different schemes 
described above. They are plotted in Fig. [21 For a total mass above ~ 0.2 — 0.3 eV the two schemes are similar and 
correspond to a degenerate scenario where each mass is ^ niiy/S. However, for smaller masses the number of neutrino 
states with relevant masses is 2 (1) in the inverted (normal) hierarchy. 

The effect of neutrino masses on cosmological observables has been usually considered equivalent for fixed ^ to^ 
(or n^h?'). However, many papers noted in the past that this is not the case and could potentially lead to differences, 
i.e. the neutrino mass spectrum should be incorporated if the sensitivity to neutrino masses is good enough (see, for 
instance the comments in 0,^^13). As an example, we note that in the mid-1990s it was shown that for CHDM 
models with the same total neutrino mass (of order some eVs), those with two degenerate massive neutrinos fitted 
better the data than those with only one (see e.g. [2^). 

Fixed the total neutrino mass, a different distribution among the 3 states {mi, TO2, tis) causes a slight modification 
of the transit from a rclativistic to a non-relativistic behavior. This can be seen in Fig. |31 where the evolution of 
the neutrino energy density is plotted for several cases with the same total neutrino mass, equally shared by 1,2 or 
3 neutrino states, as well as the realistic NH and IH schemes (taking the best-fit values of Am^). Therefore, the 
evolution of background quantities is not completely independent of the mass splitting. However, the main difference 
appears at the level of perturbations. Indeed, in the case of non-degenerate massive neutrinos, various free-streaming 
scales are imprinted in the matter power spectrum P(k). This is illustrated in Fig. ^ where we compare P(k) in the 
same cases as in Fig. O These results were obtained with our modified version of the public code CMBFAST jSO] (see 
section Ivl for details). 

We have recently summarized the effects of massive neutrinos on cosmological observables in [T^ . Here we simply 
remind that only neutrinos with masses close to the recombination temperature (Tdec 0.3 eV) leave an imprint on 
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FIG. 3: Evolution of the total neutrino energy density as a function of the scale factor of the Universe for models where the 
same ^mi/ (0.12 eV) is distributed differently. Each line corresponds to the energy density of 4 different cases (only 1 or 2 
massive states, Normal and Inverted Hierarchy) normalized to the case with 3 massive states with mass mo = ^mi//3. 
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FIG. 4: Comparison of the matter power spectrum obtained for various models where the same (0.12 eV) is distributed 

differently. The four lines correspond to the cases with 1 or 2 massive states, Normal and Inverted Hierarchy, divided each time 
by that with 3 massive states of equal mass mo = ^m,y/3. Differences in the various individual masses and free-streaming 
scales affect the position and amplitude of the break in the power spectrum. 



the CMB angular spectra, while neutrinos with smaller masses have almost the same effect as massless neutrinos. On 
the other hand, the dominant effect is the one induced by free-streaming on the matter power spectrum. Therefore, 
the usual strategy is to combine CMB and LSS measurements, where the former roughly fix most of the cosmological 
parameters, while the latter is more sensitive to neutrino masses. 



III. FUTURE CMB AND LSS DATA 



In this section wc briefly describe the experimental projects, planned or in development, that will provide data on 
the CMB anisotropy spectrum or on the distribution of LSS. 
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A. CMB experiments 

The quality of the first-year data from the Wilkinson Microwave Anisotropy Probe (WMAP) [ pi, co mplemented by 
the results of other experiments at smaller angular scales such as ACB AR, CBI or VS A |^ l32l l33l| , has shown the 
importance of CMB data as a probe of cosmological parameters. The CMB experiments measure the temperature 
fluctuations in the sky that can be expanded in spherical harmonics, 

AT 

— (0,0) =^az™rz™(0,<^). (3) 

Lni 

If the underlying perturbations are Gaussian, all information is encoded in the angular power spectrum C'l = (|a;mp). 
In addition the CMB experiments can be sensitive to polarization anisotropics, that are expressed in terms of the 
angular spectra of the E and B modes of polarization, as well as the temperature polarization cross-correlation (TE) 
spectrum. 

After WMAP, the next satellite mission will be PLANCK^, to be launched in 2007, whose experimental parameters 
are listed in TableQ] After a couple of years, it will provide CMB data more precise than that of WMAP, in particular 
concerning polarization. We also consider the CMBpol or Inflation Probe mission concept, presented in the framework 
of NASA's Beyond Einstein Program^. This experiment would have better sensitivity than the limit imposed by cosmic 
variance (up to / ~ 2300 for E-polarization, even beyond for temperature). 

In parallel to the satellite missions, there will be ground-based experiments that will measure the CMB at smaller 
angular scales with signiflcantly smaller sky coverage but good sensitivities, such as SPTpoP (in construction). ACT^ 
(funded in January 2004), or QUaD [s^l (in construction). As an example, we consider SPTpol with the characteristics 
listed in Table H] 

The observed power spectrum can be decomposed into primary anisotropics, gravitational Icnsing distortions, and 
foreground contamination. The central frequencies of CMB detectors are usually chosen in order to minimize the 
foreground contribution. In addition, by combining various frequencies, future experiments will have the power to 
separate efficiently the CMB blackbody from the various foregrounds contributions, even on small angular scales where 
the latter start to be significant. It is possible to build models for the foregrounds and to predict their impact on 
parameter extraction [34.l35l l3^: this approach is rather model-dependent, since the level of many foreground signals 
has not yet been measured experimentally. Here, we will not enter into such details. When dealing with PLANCK, we 
will employ only three frequency channels from the high frequency instrument (HFI), making the (usual) simplifying 
assumption that other channels will be used for measuring the various foregrounds, and for cleaning accurately the 
primary signal. We will do similar assumptions for SPTpol and CMBpol. We will also speculate on the results of 
an "ideal CMB experiment" limited only by cosmic variance. Then, we will limit ourselves to ^max ~ 2500 both 
for temperature and polarization, which assumes an efficient method for foreground subtraction - in particular of 
point-like sources and dust - but remains realistic (as indicated by Fig. 7 in pSjl. For the two satellite experiments, 
we assume a sky coverage of /sky = 0.65, which represents a conservative estimate of the data fraction that will be 
included in the analysis in order to avoid galactic foregrounds. For the "ideal CMB experiment" , we adopt the more 
optimistic value /sky = 1, assuming that all galactic foregrounds can be subtracted (see e.g. the component separation 
method described in ^7\). 

The issue of gravitational lensing distortion is subtle and potentially very interesting. Since lensing is induced 
by large scale structure, mainly on linear scales, this effect can be accurately predicted for a given matter power 
spectrum. Therefore, if the gravitational distortion of the CMB maps could be measured directly, there would be an 
opportunity to extract the matter power spectrum (and the underlying cosmological parameters) independently from 
redshift surveys. A way of doing this is described in [s^, E3 j and has been already applied to future neutrino 
mass extraction by |l7l |. Here, we do not incorporate this method, and assume that the matter power spectrum is 
measured only with redshift surveys, leaving a combined analysis for the future. Therefore, throughout the analysis, 
we will employ the unlensed CMB power spectra^. For the T, E and TE modes, lensing distortions are subdominant. 
In contrast, for the B mode, lensing is expected to dominate over the primary anisotropics at least on small angular 
scales. The angle above which lensing is subdominant crucially depends on the tensor-to-scalar ratio, an inffationary 



^ http : //www. rssd . esa. int /index .php?project=PLANCK 

^ http : //universe . gsf c . nasa . gov/progr a m/inflation . html ^ 

^ http://astro.uchicago.edu/spt/ 

* http : / /www. hep .upenn. edu/~angelica / act/act .html] 

Note that including lensing corrections is technically easy with CMBFAST. However, this would introduce some c orrelations among 
different modes and scales that would artificially lower the predicted errors on each cosmological parameters llTl l38l . 
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TABLE I: Experimental parameters of CMB projects: here Oi, measures the width of the beam, At,p are the sensitivities per 
pixel in ^K, v is the center frequency of the channels in GHz and /sky the observed fraction of the sky. For the PLANCK 100 
GHz channel, the value of Ap takes into account the recent design with eight polarized bolometers. 

parameter which order of magnitude is still unknown. So, we follow a conservative approach and not take the B mode 
into account. This amounts in assuming that the gravitational wave background generated by inflation is small, so 
that the B mode gives no information on primary anisotropics. 

B. Galaxy surveys 

The existing data on the distribution of galaxies at large scales come from several galaxy surveys, of which the 
completed 2dF survey® and the ongoing Sloan Digital Sky Survey*" (SDSS) are the largest. SDSS will complete 
its measurements in 2005. The matter power spectrum P{k) can be reconstructed from the data, which gives an 
opportunity to test the free-streaming effect of massive neutrinos. However, the linear power spectrum is found 
modulo a biasing factor 6^, which reflects the discrepancy between the total matter fluctuations in the Universe, and 
those actually seen by the instruments. Here we assume that the bias parameter h is independent of the scale k. 

An important point concerning LSS data is the non-linear clustering of the smallest scales. The usual approach 
is to discard any information above an effective cut-off wavenumber fcmax, while considering results at lower fc's as a 
direct estimate of the linear power spectrum. The cut-off value must be chosen with care: if fcmax is too small, we 
can lose a lot of information, especially concerning the neutrino free-streaming scale. If fcmax is too large, we can 
underestimate the error on cosmological parameters, first by neglecting any theoretical uncertainty in the quasi-linear 
corrections that could be applied to the spectrum, and second by ignoring the non-gaussianity induced by non-linear 
evolution [4ll |. 

Apart from fcmax, the important parameter characterizing the galaxy survey is its effective volume in k space, 
defined in |4ll |. If the number density of objects in the survey n(r) is roughly constant over the survey volume, and 
if the observed power spectrum P{k) is bigger than 1/n over the scales of interest (i.e., from the turn-over scale in 
P{k) up to fcmax), the effective volume is equal to the actual volume of the survey. This is a reasonable approximation 
for all the examples that we will consider here. For instance, the SDSS the Bright Red Galaxy (BRG) survey has an 
effective volume of roug hly Fcff ~ 1 (Gpc//i)3 [l| (which comes from a sky coverage /sky = 0.25 and a radial length 
of 1 Gpc /i-i). 

Beyond SDSS, plans for larger surveys are under discussion. For instance, we can mention the Large Synoptic Survey 
Telescope* (LSST), which in the future could cover the entire sky and at the same time be capable of measuring fainter 
objects LSST is designed mainly for weak lensing observations. In order to map the total matter distribution up 
to half the age of the Universe (i.e., up to a redshift z ~ 0.8 or a radial length I ^ 2.3 Gpc//i) in a solid angle 30,000 
deg^ (/sky ^ 0.75), it could measure 2 x 10* redshifts up to z = 1.5. Inspired roughly by these numbers, at the end 
of this analysis, we will speculate on the possibility to measure the power spectrum in a effective volume as large as 
Kff = (47r/3)/sky^' ~ 40 (Gpc//i)3. 

The mechanism of structure formation affects larger wavelengths at later times. So, in order to measure the linear 
power spectrum on small scales, it would be very useful to build high redshift galaxy surveys. This is one of the main 
goals of the Kilo- Aperture Optical Spectrograph (KAOS) proposal^. KAOS could build two catalogs centered around 
redshifts z = 1 and z = 3, corresponding roughly to fcmax ~ 0.2 h Mpc""'^ and fcmax ^ 0.48 h Mpc~^ respectively, 
instead of fcmax ~ 0.1 Mpc^^ today (conservative values). In both catalogs, the number density would be such that 



http : //www . mso . anu . edu . au/2dFGRS/ | 


http : //wMW . sdss . org 




http : //www . Isst . org 


http : //www . noao . edu/kaos | 
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1/n ~ -P(fcmax), and the effective volume of the two samples close to Kff ~ 0.5 (Gpc//i)'^ and VeS ~ 0.6(Gpc//i)^ 
respectively This experiment is designed mainly for measuring the scale of baryonic oscillations, in order to 
constrain dark energy. However, we will see that it would be also appropriate for improving constraints on the 
neutrino masses. 

IV. FORECAST OF FUTURE BOUNDS: FISHER MATRIX ANALYSIS 

Since the characteristics of future CMB experiments and galaxy surveys are already known with some precision, it 
is possible to assume a "fiducial" model, i.e., a cosmological model that would yield the best fit to the future data, 
and employ the Fisher matrix method to forecast the error with which each parameter will be extracted. This method 
has been widely used for many cosmological parameters, some of them related to neutrinos. For instance, we can 
mention that forecast analyses based on the Fisher matrix have shown that with future data there will be a potential 
sensitivity to an effective number of neutrinos of the order AA^off ~ 0.2 ^3,0,0], a value that is complementary to 
and will eventually improve the accuracy of primordial nucleosynthesis results (see e.g. |46l l47| 'l. 

Starting with a set of parameters Xi describing the fiducial model, one can compute the angular power spectra of 
CMB temperature and polarization anisotropies , where X = T,E,TE. Simultaneously, one can derive the linear 
power spectrum of matter fluctuations P{k), expanded in Fourier space. The error Sxi on each parameter can be 
calculated from the reduced (dimensionless) Fisher matrix Fij , which has two terms. The first one accounts for the 
CMB experiment and is computed according to ref. ji^ 

^«(CMB)='|:i:gco,-(C,-,C,>')ig^, (4) 

1=2 X.Y ■> 

where Gav{C^ ,Cj) is the covariance matrix of the estimators of the corresponding CMB spectrum. For instance, 
the TT element is given by 



Cov(Q^,Cr) = 



[21 + l)/sky 



H 2 

(5) 



ch. 



Here, the first term arises from cosmic variance, while the second is a function of the experimental parameters summed 
over the channels: = exp(— ^(Z + l)0^/81n2) is the beam window function (assumed to be Gaussian), 9b is the 
FWHM of the beam and ujt = (^bAr)~^ is the inverse square of the detector noise level (At is the sensitivity per 
pixel, and the solid angle per pixel can be approximated by 9^). For the experiments that we consider here, all these 
numbers can be found in Tabled The other terms of the covariance matrix can be found, for instance, in [l3j. 

The second term of the reduced Fisher matrix accounts for the galaxy survey data and is calculated following 
Tegmark [ill, 

f«(LSS) . 2. ^■°^o,..Wai„P...(.) 

Jo (J hi Xi o In Xj 

Here w{k) = Vcti/{2'K/kf is the weight function of the galaxy survey and we have approximated the lower limit of the 
integral fcmin — 0. We defined Pobs(fc) = b^P{k), and fcmax is the maximal wave number on which linear predictions 
are reliable. This expression is only an approximation, since in addition to non-linear clustering it ignores edge effects 
and redshift space distortions. 

Inverting the total Fisher matrix, one obtains an estimate of the l-cr error on each parameter, assuming that all 
other parameters are unknown 

^ = iF-')\l'- (7) 

Xi 

It is also useful to compute the eigenvectors of the reduced Fisher matrix (i.e., the axes of the likelihood ellipsoid in 
the space of relative errors). The error on each eigenvector is given by the inverse square root of the corresponding 
eigenvalue. The eigenvectors with large errors indicate directions of parameter degeneracy; those with the smallest 
errors are the best constrained combinations of parameters. 



The characteristics of KAOS are taken from the "Purple Boole" available on-line at |http : //www . noao ■ edu/kaos | 
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FIG. 5: Predicted 2a error on the total neutrino mass M = ^ as a function of M in the fiducial model, using PLANCK 
and SDSS (limited to fcmax = 0.15 h Mpc^^). The left plot was obtained with the preferred experimental value of Amatmi a^nd 
the right plot with the current 3a" upper bound. In each case, we show the results assuming either NH or IH. 

V. RESULTS 

We have computed the total Fisher matrix from Eqs.EJandEl using various experimental specifications. Throughout 
the analysis, our fiducial model is the concordance "flat ACDM" scenario, with parameters close to the current best- 
fit values and with additional neutrino masses. The nine free parameters with respect to which derivatives are 
computed arc: flmh^ (matter density, including baryons, cold dark matter and neutrinos), fli,h^ (baryon density), 
r^A (cosmological constant), CJqq (amplitude of temperature spectrum at multipole 200), Ug (scalar tilt), r (optical 
depth to reionization), i/hp (fraction of baryonic mass in the form of Helium), M = '^m^ (total neutrino mass) and 
b (unknown bias of the LSS data). The fiducial value of b is irrelevant by construction, and we will try various values 
of M, distributed following the NH or IH scheme. Other fiducial values read: 

{n^h^,nbh^ ^A, Clm^n,,T, yne) = (0.143, 0.023, 0.70, 0.85, 0.96, 0.11, 0.24). 

All derivatives are computed at zero spatial curvature (by varying h appropriately). Note that we use double-sided 
derivatives with step 10% for M, 50 % for yne, 5% for all other parameters. We checked carefully that these steps 
are sufficient in order to avoid possible numerical errors caused by the limited precision of the Boltzmann code - in 
our case, version 4.5.1 of CMBFAST [^OIj with option "high precision". We also checked that with twice larger steps, 
the results change only by a negligible amount. These conditions were not a priori obvious for the smallest neutrino 
masses studied here, but we increased the precision of the neutrino sector in CMBFAST accordingly. Actually, in order 
to study three neutrino species with different masses, we performed significant modifications throughout CMBFAST. 
For each mass eigenstate, we integrate some independent background and perturbation equations, decomposed in 
15 momentum values, up to multipole / = 7. Finally, we include the small distortions in the neutrino phase-space 
distributions caused by non-instantaneous decoupling from the electromagnetic plasma (with QED corrections at 
finite temperature) |2g, but these last effects are almost negligible in practice. 

A. PLANCK+SDSS 

We first derive the precision with which the combined PLANCK and SDSS data will constrain the total neutrino 
mass in a near future. Experimental specifications for these experiments are given in the previous section, and we 
choose to limit SDSS data to the scale fcmax = 0.15 h Mpc~^ where non-linear effects are still small. Fig. |51 shows 
the predicted 2a error on M for various fiducial models, assuming different values of M, the two possible schemes for 
the mass splitting (either NH or IH), and two different values of Am^^^. The solar mass scale Am^^^ is essentially 
irrelevant in this analysis, and is kept fixed to the current preferred value of 6.9 x 10^^ cV^. The possible values of 
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TABLE II: Absolute errors at the X-a level for various cosmological models, using PLANCK+SDSS (fcmax ~ 0.15 /i Mpc ^). 
The first line shows our simplest flat ACDM model, described by 9 free parameters with fiducial values C200 = 0.85, Us — 0.96, 
T = 0.11, = 0.70, n^h^ = 0.143, nth^ = 0.023, M = 0.3 eV (normal hierarchy), Yne = 0.24. The value chosen for 
b^P{ko = 0.1 h Mpc~^) is irrelevant. The next lines have one additional parameter X: an effective number of neutrinos A*'^ 
parametrizing the abundance of extra relativistic relics, with fiducial value 0; a free spatial curvature parametrized by r2k with 
fiducial value 0; a free time-independent equation of state for dark energy parametrized by w with fiducial value —1; a free 
scalar tilt running parametrized by a = dus/dlnk with fiducial value 0. 



M are of course bounded from below: the minimal value corresponds to the limit in which the lightest neutrino mass 
goes to zero, in one of the two NH or IH schemes. 

Let us first concentrate on the case in which Am^^j^^ has its current preferred value of 2.6 x 10"'^ eV^ (left plot). The 
minimal value of M in the NH (resp. IH) case is approximately 0.06 eV (resp. 0.10 eV). However, the 2(t detection 
threshold, defined by M = 2a{M), is around 0.21 eV. We conclude that PLANCK+SDSS will probe mainly the region 
were the three neutrinos are quasi-degenerate in mass, with no possibility to distinguish between the two cases. In 
absence of clear detection, the 2a upper bound will be of order 0.2 eV, corresponding to individual masses (0.08, 0.06, 
0.06) eV assuming NH, or (0.073, 0.073, 0.053) eV assuming IH. As expected, we find that the 2a detection threshold 
is still 0.21 eV when the calculations are performed with a larger value Ato^j,jj = 3.7 x 10~^ eV^ (the 3a upper bound 
in Eq. as shown in the right plot of Fig. |S1 

It is interesting to study whether this precision is limited mainly by a degeneracy between Af and some combination 
of other cosmological parameters, or simply by the experimental sensitivity to the individual effect of M . In the first 
case, the results could be improved by including priors from other types of experiments on the cosmological parameters; 
in the second case, one would have to wait for a new generation of CMB and/or LSS experiments. In order to address 
this point, we computed the eigenvectors and eigenvalues of the reduced Fisher matrix. It turns out that for all 
our fiducial models, one of the unit eigenvectors points precisely in the direction of M , with coefficient very close to 
one in this direction (and, of course, the corresponding eigenvalue matches the error previously obtained for M). We 
conclude that M is not affected by a parameter degeneracy, and that independent measurements of other cosmological 
parameters would not help very much in constraining neutrino masses. Note that this is not yet the case for current 
cosmological bounds on neutrino masses, where the addition of priors on parameters such as the Hubble constant or 
r^A leads to more stringent bounds (see e.g. [l^). 

The absence of large parameter degeneracies applies to our reference model with nine free parameters. It may not 
necessarily be true in the presence of extra parameters describing deviations from the concordance ACDM model. 
In order to illustrate this point and to test the robustness of our conclusions, we have calculated the error on each 
parameter for several extended cosmological scenarios, with extra relativistic degrees of freedom, spatial curvature, 
dark energy with varying density but constant equation of state, or a primordial spectrum with running tilt (see Table 
HTj) . The neutrino mass bound is found to be quite robust in all these cases, which proves that in none of these models 
the effect of M can be mimicked by some other parameter combination. 

It is also interesting to study the relative impact of CMB temperature, CMB polarization and LSS data on the 
measurement of M. We show in table IIIII the error on each parameter for SDSS alone, PLANCK alone (with or 
without polarization), and various combinations of CMB and LSS data, with an explicit dependence on the value of 
^max- The complementarity of PLANCK and SDSS clearly appears. While PLANCK alone would achieve only a la 
detection of M = 0.3 eV and SDSS alone would not detect it at all, the combined data would probe this value at the 
3a level. One can check from Table Hill that PLANCK data on polarization lowers the error on M by approximately 
30%. By diagonalizing the "PLANCK (no pol.)-|-SDSS" Fisher matrix, we checked that without polarization there 
would be a significant degeneracy between neutrino mass and optical depth to reionization. Indeed, while reionization 
lowers the CMB temperature spectrum keeping the matter power spectrum unchanged, the effect of neutrino free- 
streaming is opposite in first approximation (at least on small scales). So, polarization measurements are indirectly 
a key ingredient for neutrino mass determination. 
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fcmax (fe/Mpc) 


In C200 


ris 


T 




n /)^ 




M (eV) 




ln[&2p(fco)] 


SDSS alone 


0.10 


- 


0.6 


- 


0.8 


0.5 


0.1 


7.0 


- 


0.3 




0.15 


- 


0.5 




0.09 


0.4 


0.08 


1.5 


- 


0.06 




0.20 


- 


0.1 


- 


0.05 


0.09 


0.02 


0.5 


- 


0.01 


PLANCK (no pol.) 


- 


0.005 


0.02 


0.10 


0.05 


0.006 


0.0006 


0.42 


0.03 


- 


PLANCK (no pol.) + SDSS 


0.10 


0.005 


0.02 


0.08 


0.02 


0.002 


0.0004 


0.24 


0.02 


0.015 




0.15 


0.005 


0.02 


0.08 


0.01 


0.001 


0.0003 


0.15 


0.02 


0.008 




0.20 


0.005 


0.01 


0.07 


0.006 


0.0009 


0.0003 


0.13 


0.02 


0.005 


PLANCK (all) 


- 


0.005 


0.008 


0.005 


0.04 


0.004 


0.0003 


0.30 


0.01 


- 


PLANCK (all) + SDSS 


0.10 


0.005 


0.007 


0.005 


0.02 


0.002 


0.0002 


0.19 


0.01 


0.012 




0.15 


0.005 


0.007 


0.005 


0.01 


0.001 


0.0002 


0.11 


0.01 


0.007 




0.20 


0.004 


0.007 


0.005 


0.006 


0.0008 


0.0002 


0.08 


0.01 


0.005 


CMBpol 




0.003 


0.003 


0.003 


0.006 


0.0006 


0.00008 


0.07 


0.004 




CMBpol + SDSS 


0.10 


0.003 


0.003 


0.003 


0.006 


0.0006 


0.00008 


0.07 


0.004 


0.011 




0.15 


0.003 


0.003 


0.003 


0.005 


0.0006 


0.00007 


0.06 


0.004 


0.006 




0.20 


0.003 


0.003 


0.003 


0.004 


0.0005 


0.00007 


0.05 


0.004 


0.004 



TABLE IIL Absolute errors at the 1-cr level, for various experiments and the same ACDM model as in table [H] (with 9 free 
parameters). In particular, the fiducial value of the total neutrino mass is still M — 0.3 eV. When using SDSS, we show the 
results for three choices of fcmax, the maximal wavenumber on which the data are compared with linear theory predictions: 
fcmax = 0.10 ft Mpc~^ (conservative), 0.15 ft Mpc~^ (reasonable), or 0.20ft Mpc~^ (optimistic). 





SDSS 


SDSS+KAOS 


"hypothetical LSS" 


PLANCK 


0.21 


0.16 


0.11 


CMBpol 


0.13 


0.10 


0.09 


"ideal CMB" 


0.10 


0.09 


0.08 



TABLE IV: 2-CT detection threshold (in eV) for various combinations of CMB and LSS experiments (assuming the normal 
hierarchy scenario). The "ideal CMB" experiment is limited only by cosmic variance up to multipole I = 2500 and covers 100% 
of the sky. The "hypothetical LSS" survey has a volume VcH — 40 (Gpc/h)^ and probes the linear spectrum up to fcmax ~ 0.15 ft 
Mpc"^ (that would be the case of a large galaxy survey covering 75% of the sky up to 2 = 0.8). 



B. Post-PL ANCK experiments 

Here we consider whether future CMB and LSS experiments will reach a better sensitivity on the neutrino mass, in 
particular at the level of the small values of M expected for the hierarchical normal and inverted schemes. Sensitivities 
significantly better than 0.1 eV would mean approaching the absolute minimum of M in the NH case or even ruling 
out the IH scenario. 

In the previous section, we mentioned a few CMB missions that have been proposed so far in complement to 
PLANCK. We will study the impact of a few of them, and of an "ideal CMB experiment" that would be limited only 
by cosmic variance up to I ~ 2500 (both for temperature and polarization). The main difficulty for reaching this goal 
would be to subtract accurately small-scale foregrounds, and in particular point- like sources, but even with current 
technology such an ideal experiment is conceivable. On the other hand, it is difficult to specify the characteristics of an 
ideal LSS experiment, since it will be limited by technological improvements in instrumentation and data processing. 
Therefore, we will keep in the analysis a free parameter VcS describing the effective volume of an ideal volume-limited 
survey. 

We show in Fig. |H|the predicted 2a error in four cases corresponding to SPTpol (upper left), PLANCK (upper 
right), CMBpol (lower left), and our ideal CMB experiment (lower right). The value of 2(t (in eV) is shown with 
grey levels, as a function of fcmax (horizontal axis) and Vcs (vertical axis) in units of T4ff (SDSS) — 1 (Gpc/h)^. The 
total mass has been fixed to M = 0.11 eV, distributed according to the NH scheme. We learned from the previous 
subsection that for higher values of M, the error could be smaller (at most by a factor 2). However, we are now 
interested in the range 0.05 eV < M < 0.2 eV, since larger values should be detected by PLANCK-I-SDSS, and 
smaller values are excluded by oscillation experiments. In this range, on can safely interpolate the results obtained at 
M = 0.11 eV. In particular, our results for a cosmic- variance limited CMB experiment are in reasonable agreement 
with those of [l5l |. 
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FIG. 6: The grey regions are the 2a expected errors on "Y^rriu (eV) for a fiducial value of 0.11 eV, as a function of the 
parameters of the galaxy survey, where each panel corresponds to a specific CMB experiment. The vertical lines indicate the 
cut-off wavenumber fcmax for the linear matter power spectrum at the conservative (optimistic) value 0.15(0.2) h Mpc~^. The 
thin contours shown are (from bottom to top) for 0.3 and 0.2 eV, while the thick contours correspond to the minimum values 
of ^ nil, in the IH (lower lines) and NH (upper lines) schemes, assuming the best-fit (solid) or the 3(t upper bound (dashed) 
value of Amatm- 

For SDSS (or for any survey with z < 1) we expect the relevant value of /cmax to be around 0.15 Mpc~^. 
However, depending on the overall amplitude of the matter power spectrum (often parametrized by erg, and still 
poorly constrained) and on future improvements in our understanding of non- linear corrections, this value might 
appear to be either too optimistic or too pessimistic: this is the reason why it is interesting to leave it as a free 
parameter. 

One can see that replacing PLANCK by CMBpol would lead to a better sensitivity to the neutrino mass, with a 
2(7 detection threshold at 0.13 eV instead of 0.21 eV. The expected errors for CMBpol, with and without SDSS data, 
can be found in Table Hill Adding to SDSS the two KAOS surveys (centered around z = 1 and z = 3) would also 
lead to some improvement. For Planck+SDSS-|-KAOS we get a 2a detection threshold of M ~ 0.16 eV, while for 
CMBpol+SDSS+KAOS one could reach M ~ 0.10 eV. These results are summarized in Table Hvl 

There is still room for improvement beyond this set of experiments. In order to make a precise statement on the 
conclusions that could be drawn on the long term, we keep the "ideal CMB experiment" characteristics and fix Ves 
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Am\tm= 2-6 X 1 0"^ eV^ (best fit) 




0.1 0.2 

fiducial M (eV) 



0.7 
0.6 
0.5 
0.4 
0.3 
0.2 
0.1 



Am\^^= 3.7 X 1 0"^ eV^ (+3o bound) 




0.05 





NH 






IH 

/// - 















0.1 0.2 

fiducial M (eV) 



0.5 



FIG. 7: Predicted 2a error on ^ as a function of ^ in the fiducial model, using an ideal CMB experiment (limited 
only by cosmic variance up to Z = 2500, both for temperature and polarization) and a redshift survey covering 75% of the 
sky up to z ~ 0.8 (Kff = 40 (Gpc//i)^), still limited to fc^ax = 0.15 ft Mpc~\ The left plot was obtained with the preferred 
experimental value of Amatmi and the right plot with the current 3ct upper bound. In each case, we show the results assuming 
either NH or IH. 



to 40 (Gpc//i)3 (in section ITIll wc argued that this could hopefully represent the volume of a survey comparable to 
the LSST project), while keeping fcmax = 0-15 ft. Mpc^^. In Fig.0 we plot the corresponding results in the same way 
as we did for PLANCK+SDSS. Assuming the IH scenario, we see that any value of the mass could be detected at the 
2ct level. Assuming NH, this is only true at the 1 or l.Scr level, depending on the value of Am^^jj^. The 2ct detection 
threshold is at 0.08 eV. 

Our results show, for the first time, that if the available cosmological data are precise enough, the expected errors 
on the neutrino masses depend not only on the sum of neutrino masses, but also on what is assumed for the mass 
splitting between the neutrino states. As can be seen from Figs. and [7| the sensitivity on M will be slightly better 
in the NH case in the mass region close to the minimum value of the IH scheme. These small differences arise from 
the changes in the free-streaming effect that we have described in section^ and obviously disappear for a total mass 
in the quasi-degenerate region (above 0.2 eV or so). 

In any case, the main contribution of cosmology to the possible discrimination between the neutrino mass schemes 
will still be the possibility of ruling out the case in which the masses are quasi-degenerate. Even in our most optimistic 
forecast (Fig. [Tj), if the preferred value of M turns out to be smaller than 0.1 eV, the error bar will still be too large 
in order to safely rule out the IH case. We also performed an extended analysis in which, instead of assuming either 
normal or inverted hierarchy, we introduced a tenth free parameter accounting for a continuous interpolation of the 
mass spectrum between the two scenarios, for fixed M . By computing the error on this parameter, we obtained a 
confirmation that the NH and IH scenarios cannot be discriminated directly from the data. However, any analysis of 
future, very precise cosmological data must take into account the texture of neutrino masses in order to translate the 
corresponding positive signal (or bound) into M. 



VI. CONCLUSIONS 



In this paper we have analyzed the sensitivities of future CMB and LSS data to the absolute scale of neutrino masses, 
taking into account realistic experimental sensitivities and extending the results of previous works 0, 0, 0, . 

We have considered the values of neutrino masses distributed according to the presently favored three neutrino 
mass schemes, that follow either a normal or an inverted hierarchy. As discussed in section^ a different distribution 
of the same total neutrino mass leads to small changes in the cosmological evolution of neutrinos, and in particular 
in the free-streaming scales (qualitatively discussed, for instance, in 0,0)- These changes disappear when the total 
neutrino mass enters the quasi-degenerate region. 

We used the Fisher matrix method to forecast the errors on cosmological parameters that can be extracted from 
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future CMB experiment and redshift survey data, assuming a fiducial 9-dimensional cosmological model close to the 
currently favored ACDM model. Our theoretical CMB and matter power spectra were generated with the standard 
Boltzmann code CMBFAST, modified in order to include three neutrino states with different masses. 

In particular, for the case of PLANCK and SDSS we found good agreement with the results of with a 2tT-error 
on the total neutrino mass of 0.2 eV that will allow us to probe only the quasi-degenerate neutrino mass region. 
Better sensitivity will be achieved with the combination CMBpol and SDSS, for which we foimd 0.12 eV, close to the 
minimum value of the total neutrino mass in the inverted hierarchy case. These results correspond to a conservative 
value of fcinax = 0.15 ft. Mpc~^, the maximal wavenumber on which the LSS data are compared with the predictions of 
linear theory. We also tested how the errors change when including additional cosmological parameters to our fiducial 
model. In general, we found that the errors on the neutrino masses are not modified in a significant way. 

Our results show that the approach where CMB experiments are only limited by cosmic variance (as in 
probably too simplistic. However, if a future CMB experiment is capable of getting close to such an ideal limit, then 
the combination with data from galaxy redshift surveys larger than SDSS would lead to errors on the total neutrino 
mass comparable to the minimum values of the hierarchical scenarios. In such a case, we have shown that there exist 
slight differences in the expected errors between the two hierarchical neutrino schemes for the same total neutrino 
mass. 

In conclusion, we consider that cosmological data can provide valuable information on the absolute scale of neutrino 
masses, that nicely complements the present and future projects of beta decay and neutrinoless double beta decay 
experiments. This conclusion is reinforced when one takes into account other cosmological probes of neutrino masses, 
complementary to the approach of the present paper. We can cite, for instance, studies of the distribution of matter in 
the Universe through the distortions of CMB maps caused by gravitational lensing (measured from non-gaussianities in 
the CMB maps) ^3 ^^'^ the weak gravitational lensing of background galaxies by intervening matter [3 Isol [5ll | . 

It is interesting to note that any information on the absolute neutrino masses from cosmology will be interesting not 
only for theoretical neutrino models, but also for connected baryogenesis scenarios which occur through a leptogenesis 
process (see e.g. [H^IH^IIH^ V 
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